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The long-standing deficit of measured versus predicted solar neutrino fluxes is re- 
examined in light of possible astrophysical solutions. In the last decade, solar 
neutrino flux and helioseismic measurements have greatly strengthened the case 
for non-astrophysical solutions. But some model-independent tests remain open. 



The solar neutrino problem has nagged physicists for over 30 years and, 
for most of us, has a natural solution in neutrino oscillations, either in vacuo 
or matter-enhanced (MSW effect). This modification of neutrino properties 
(with m v <, 10 -2 eV and mixing ~ 0.001-1) requires an extension of the Stan- 
dard Model not detectable in accelerators and consistent with many models 
of unification, u But confidence in such a solution is based on a prior exclusion 
of astrophysical or nuclear physics explanations; only in the past decade has 
such an outcome become strongly credible. 

Predictions of the solar neutrino flux <fi are outputs of a solar model, which 
in turn are special cases of hydrogen-burning (4-H — > 4 He) main sequence 
stellar models (almost all pp chain, with small CNO contribution). These 
predictions are usually quoted from a particular model (here the Bahcall- 
Pinsonneault 199J*Lor BP98 model), and detailed models agree if the same 
inputs are used. Era But a more generic approach is attractive if it can free 
us from a specific model with fixed parameters. Greater generality is all the 
more important if it reveals basic properties of the Sun and <f> that depend 
only on simple properties. Here I outline the results of such an approach. nUu 
As far as observations now go, it confirms the results of detailed solar models 
but also specifies crucial solar observations that remain to be filled in. 

1 Properties of Solar and Stellar Models 

Simple solar models orient us with gross structure: the core {t/Rq 0.3, 
where nuclear fusion generates the luminosity L©), the radiative zone (0.3 
& r/R® < 0.71) and the convective zone (CZ, r/R @ > 0.71) (Figure [l]). Less 
massive stars (M < Mq) have, according to stellar models, even deeper con- 
vective zones; while for M > Mq, the outer convective zone disappears, and 
a convective inner core appears as the central temperature gradient surpasses 
a critical value. The Sun might have had a convective inner core, drastically 



HepOO-7: formatted in World Scientific on February 5, 2008 



Figure 1. Schematic cutaway of the Sun in a generic standard solar model. 



changing the <f> predictions, but there is now decisive evidence against core 
convection (see Section || below). A higher central temperature would also 
have led to CNO-dominance and higher neutrino fluxes. El 

A solar model is standard (an SSM) if the model contains all the physics of 
matter, gravitation, and nuclear fusion needed to obtain a star, but nothing 
more. Conceptually, stellar structure and evolution divide into three levels 
by time scales. For the Sun, chemical evolution needs ~ 10 Gyr; thermal 
equilibrium, about 10 Myr (Kelvin- Hclmholtz time); and hydrostatic equilib- 
rium, about 5 minutes. The hydrodynamic time controls the helioseismic p- 
and g-modes. (Late in the Sun's evolution, the chemical time scale will be 
shortened and the hierarchy bbirred.) Structure means only the thermal and 
mechanical features of a star, vffl 

The properties of matter needed are constitutive relations giving pressure 
P, opacity k, and specific luminosity generation e as functions of density p and 
temperature T . If the thermal structure is given, then the equation of state 
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reduces to a barytrope P — P(p), and the mechanical structure alone becomes 
a closed problem characterized by the stiffness profile r = d In P/d In p. The 
special case ofr = 7 = l + l/n = constant is a polytrope of index n. All n 
polytropes have finite mass and radius, n = is the constant-density case, clu 
The initial conditions are fixed mass M and the element mass abundances 
Xi (the zero-age star assumed chemically homogeneous) . The boundary con- 
ditions are zero mass and luminosity at the center and (nearly) zero pressure 
and density at the surface. The Xi develop gradients by evolution, as nuclear 
fusion and-heavy element diffusion act, and additional helium accumulates in 
the core. HO 

Even within the SSM framework, variations are possible. The k and s 
functions and fusion cross sections must be calculated from atomic and nuclear 
physics and extrapolated into regimes not directly testable. For the dominant 
luminosity-producing ppl reactions (terminating through 3 He- 3 He fusion to 
4 He), e and the reaction rates are almost fixed by Lq. But reactions not 
strongly connected with Lq are not well-constrained: the ppl I and ppIII 
chains (terminating from 7 Be through 7 Li and 8 B to 2 • 4 He, respectively) are 
sensitive functions of the core temperature and nuclear cross sections, as well 



as mildly dependent on the core density. Lq and othe 
only weak constraints on the ppl I and ppIII rates. [ 



obal properties place 



2 Generalizing the Standard Solar Model 

Generic properties of solar structure are restricted by boundary conditions. 
The SSM can then be generalized in two ways. One is to calibrate with the 
specific model. This procedure defines a generalized SSM family, although 
not the most general. 

Its most convenient implementation is through homology or power-law 
scaling, derivable analytically or made evident by numerical solutions. Ex- 
ploration of model space by varying SSM inputs is actually a special case 
of homology, which amounts to "small" perturbations of the logarithms of 
inputs and outputs. Such "perturbative" analysis works ff¥£ T a surprisingly 
large range, so long as the power law-relations are stable. Cm 

The first signs of homological behavior iruSSMs were found in the 1000 
SSM Monte Carlo study of Bahcall and Ulrich. H Subsequent work over a wider 
model range revealed a much broader validity for homology. rThe underlying 
analytic structure was derived by Bludman and Kennedy, u Starting with 
structure alone, one keeps only dimensional and scaling behavior of macro- 
scopic variables, dropping the differential nature of the equations. Assum- 
ing multifactor power laws for the equation of state, opacity, and luminosity 
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generation, we found homological relations for the mechanical and thermal 
structure, assuming fixed powers. This requirement restricts the homology 
to the radiative and core regions, below the CZ. The dominant luminosity 
production is by ppl, carried outwards entirely by radiative diffusion. The 
constitutive relations are 

P/p = WT/ l i , K (p,T) = K (X l )p n T- s , e{p,T) = e {X)p x T v . (1) 

Expanded about the SSM, the exponents are: n = 0.43, s = 2.5, A = 1.0, v 
= 4.2. jji. re G , and e G are composition-dependent. The luminosity constraint 
reads: Era 

</>(pp) + (O.977)0(Be) + (O.751)0(B) + (O.956)0(CNO) = 6.55xl0 10 cm" 2 sec" 1 

(2) 

The boundary conditions are imposed in a way appropriate to a sin- 
gle star: Mq, L@, and Rq fixed. Homology then gives a family of possi- 
ble non-convective interiors consistent with observed outer solar features and 
parametrized by p c and T c : 

--0. 34 -0.40 ,,0.52 r 0.085 rp ^-0.13 -0.034,0.22 rO. 17 /o\ 

p c ~ e Q k p c Liq , i c ~e Q k p c . [3) 

The resulting 4> scale stably with p c and T c over a large range, giving 
the two-parameter homological mechanical/thermal variations of the SSM: 
(j>(i) ~ p^ -Tf% with (o^A) for pp, Be, and B u'b being (-0.1, -0.7), (0.7,9), 
and (0.3, 21), respectively. The highest reactions in the pp chain have the fa- 
mous extreme sensitivity to T c , while all sensitivities to p c are mild and arise 
from the small luminosity contribution made by the ppl I, ppl II, and CNO 
chains. It should be stressed that p c and T c are model outputs, like the 
These exponents reproduce the 1000-SSM Monte Carlo and clarify that the 
entire homological class of SSMs has the wrong pattern of fluxes to explain 
the observed energy dependence of (f> : lower energies are more suppressed. 
(Variation of nuclear cross sections also fail to explain the pattern.) This 
conclusion depends only on pp/-dominance in e, radiative diffusion in re, and 
the ideal gas law. 

It is instructive to compare these homology results with the approxima- 
tions often used to model stars. Perhaps the simplest (after the constant- 
density case) is the Eddington standard model, based on a constant ratio of 
radiation to matter pressure throughout the star and equivalent to an (n, 7) 
= (3, 4/3) polytrope. u This model, once its free parameters are fit, represents 
many main sequence stars not badly. Homology applied to the mechani- 
cal structure alone automatically leads to a polytrope. Q But the Bludman- 
Kennedy homology is more general than a polytrope, as it applies to both 
mechanical and thermal structure. It also scales correctly in the evolved core, 
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Figure 2. Stiffness profile T(r) in the Bahcall-Pinsonneault 1998 SSM. |'| 



where the molecular weight (i changes substantially, reflecting the chemical 
evolution that makes the present Sun differ from its zero-age incarnation. No 
polytrope fits this behavior. uB 

A more complete version of homology is possible if the differential struc- 
ture is retained and rewritten using scale-invariant homology variables. B For 
the entire mechanical, thermal, and chemical structural system, the dimen- 
sionless differential equations are not less complex than a full SSM. But if the 
structure is restricted to the mechanical alone and a barytrope P(p) assumed, 
simple dimensionless structure equations follow. The key to the mechanical 
structure turns out to be the T profile (Figure ||). Constant T gives a poly- 
trope again; in fact, Tssm — 4/3 outside the inner core, up to the CZ, where 
it rises to 5/3 (the adiabatic value). But within the core, \x rises and T drops 
towards to the center, where r^ M ~ 8/9. Described in terms of a polytrope, 
the effective index n c s = (T — rises in the core, diverging at T = 1. 
Further towards the center, n e ff rises from minus infinity to a finite negative 
value at the center. Such behavior in the inner core is not even approximately 
polytropic and explains why attempts to use polytropes to approximate the 
Sun only work (and only crudely) over the whole Sun and fail badly in the 
core. 
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Figure 3. Bouyancy frequency profile N(r) in the Bahcall-Pinsonneault 1998 SSM. 



The other approach to generalizing the SSM is to work with a few reason- 
able assumptions, following these to simple, testable predictions. Some power- 
ful results are available, although restricted to mechanical structure only, □ to 
which helioseismology is the key: adiabatic sound waves are mechanical per- 
turbations with information about sound speed, equation of state, and density 
and pressure profiles. 12.3 Taking full advantage of these results requires both p- 
(pressure) and g- (gravity) modes. Their spectra can be inverted to yield adi- 
abatic sound speed c a( j(r) and bouyancy frequency N(r) profiles, from which 
follows the complete mechanical structure. The BP98 N(r) profile is shown 
in Figure ||. 

The thermal and chemical structures cannot be directly probed by he- 
lioseismology, as their associated time scales are so long. Simple homology 
implies L ~ /i 4 M 3 k^ 1 e~° , but changes in Lq large enough to explain the solar 
v deficit are probably ruled out by paleoclimatology. Eil Direct tests of these 
aspects of the SSM require comparison to other Sun-like stars. Such stars 
vary from the Sun somewhat in mass and chemical composition and could 
lie anywhere on their respective evolutionary tracks. Comparison properties 
include luminosity, surface temperature, and photospheric radius. With accu- 
rate photometry and parallaxes, accurate luminosities and colors are achiev- 
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Figure 4. Solar sound speed profile c a( j(r) inferred from helioseismic observations 



able. Il3 The intermediary between these observations and stellar structure is 
stellar atmosphere models, which have advanced considerably in the last 30 
years. Although still oversimplified, the models are good enough for .solar- 
type stars to infer ranges for surface T, g, abundances, and turbulence. LL3 An 
exciting possibility will be opened by asteroseismology of Sun-like stars, as 
observation of stellar seismic modes ^especially g-modes) would lead to direct 
characterization of stellar interiors. lia 



3 Observational Issues 

Mq, Lq, i?0, and surface T, as well as surface and proto-solar (meteoric) 
abundances, are well measured. Helioseismic observations have directly or in- 
directly captured millions of p-mod£s, allowing an accurate inversion of c a( j(r) 
down to t/Rq = 0.05 (Figure ^). Ej A convective core (7 = 5/3) is ruled out, 
although circulation of heavy elements not affecting heat transport cannot be 
at present. Ej 

The inferred sound speed peaks off-center at r/R & ~ 0.07. Hydrostatic 
equilibrium requires dc^/dr = at the center, but an off-center peak occurs 
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where T = 1. □ This peak and dc^/dr > for r/R Q < 0.07 indicate T < 1 
there, a crucial confirmation of the core's chemically evolved state. A com- 
plete profile down to the center becomes possible with the lowest p-modes. 
But complete inversion for model-independent mechanical structure would be 
possible only if the higher g-modes were observed; analogous inversion would 
yield N(r), and together N and c a( j yielcLT and other mechanical profiles, the 
first truly independent test of the SSM. □ 

Comparing the Sun with other sun-like stars has been possible for many 
decades, albeit at poor precision. But the recent Hipparcos-Tycho star cat- 
alogs (edition I in 1997: 1.1 M stars; edition II in 2000: 2.5 M stars) have 
revolutionized astrometry, raising the accuracy .of nearby stellar parallaxes 
by up to a factor of 10 or better (1 m-arcsec). Il3 Luminosities accurate to 
< 5% for nearby stars (within 25 pc) can now be inferred. With good color 
measurements and_b£st current stellar atmosphere models, surface T's can 
be limited to 1%. E2Hiij Even more dramatic astrometric improvements could 
come from the proposed SIM and GAIA orbital systems, to be launched in 
2006 and 2009, respectively: 4 /x-arcsec parallax errors and luminosity errors 
limited only by photometry, tenths of a percent. 113 
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